Cosmic-ray and X-ray Heating of Interstellar Clouds and 

Protoplanetary Disks 



Alfred E. Glassgold 

Astronomy Department, University of California, Berkeley, CA 94720-3411, USA 

aglassgold@berkeley . edu 
Daniele Galli 

INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125, Italy 

galliOarcetri . astro . it 
and 

Marco Padovani 

Laboratoire de Radioastronomie Millimetrique, UMR 8112 du CNRS, Ecole Normale 
Superieure et Observatoire de Paris, 24 rue Lhomond, 75231 Paris cedex 05 

marco . padovani@lra . ens . f r 

Received 2012 May 14; accepted 2012 July 20 



- 2 - 



ABSTRACT 

Cosmic-ray and X-ray heating are derived from the electron energy loss cal- 
culations of Dalgarno, Yan and Liu for hydrogen-helium gas mixtures. These 
authors treated the heating from elastic scattering and collisional de-excitation 
of rotationally excited hydrogen molecules. Here we consider the heating that 
can arise from all ionization and excitation processes, with particular emphasis 
on the reactions of cosmic-ray and X-ray generated ions with the heavy neutral 
species, which we refer to as chemical heating. In molecular regions, chemical 
heating dominates and can account for 50% of the energy expended in the cre- 
ation of an ion pair. The heating per ion pair ranges in the limit of negligible 
electron fraction from ~ 4.3 eV for diffuse atomic gas, to ~ 13 eV for the moder- 
ately dense regions of molecular clouds and to ~ 18 eV for the very dense regions 
of protoplanetary disks. An important general conclusion of this study is that 
cosmic-ray and X-ray heating depends on the physical properties of the medium, 
i.e., on the molecular and electron fractions, the total density of hydrogen nuclei, 
and to a lesser extent on the temperature. It is also noted that chemical heat- 
ing, the dominant process for cosmic-ray and X-ray heating, plays a role in UV 
irradiated molecular gas. 

1. Introduction 

There has been a significant increase of interest in the cosmic-ray ionization rate in both 
local and distant neighborhoods. Much of it has been spurred by infrared absorption-line 
measurements of the Hjj" ion in diffuse and translucent interstellar clouds (e.g., Indriolo et 
al. 2007 and references therein). Indriolo & McCall (2012) have extended these measurements 



-3 - 



to fifty lines of sight dominated by translucent clouds with more than one magnitude of visual 
extinction, obtaining 21 detections. Their straightforward analysis of the observations yields 
values of the ionization rate for H 2 in the range ~ (2 — 11) x 10~ 16 s -1 , significantly larger 
than the range of values that had been in use for several decades^. This result has led to the 
reconsideration of how both electron and proton cosmic-rays propagate in the interstellar 
medium (e.g., Indriolo et al. 2009, Padovani et al. 2009, Padovani & Galli 2011, Rimmer 
et al. 2012, Everett & Zweibel 2011). In the context of the new survey, Indriolo & McCall 
al. (2012) reviewed the explanation advanced by Indriolo et al. (2009) and Padovani et 
al. (2009) that the intensity of low-energy protons is reduced inside dense molecular clouds. 
The large difference between diffuse/translucent and dense clouds is maintained in Indriolo 
& McCall (2012), but their sample shows little evidence for a dependence of ionization rate 
on column density. The value of the cosmic ray ionization rate for the local interstellar 
medium remains an open question. 

Another topic of interest relating to the interaction of cosmic-rays with the interstel- 
lar medium is their role in heating. Cosmic-rays are an efficient (often dominant) source 
of heating in various environments, from the dense gas in molecular clouds (Goldsmith Sz 
Langer 1978), both in normal and starbust galaxies (e.g., Suchkov et al. 1993), to pho- 
todissociation regions (Shaw et al. 2009), and possibly even in the primordial gas (Jasche 
et al. 2007). In dense, shielded regions like molecular cloud cores, the balance of cosmic- 
ray heating and gas cooling determines a temperature gradient decreasing inwards (Galli et 
al. 2002) that has been accurately traced by interferometric observations of molecular emis- 
sion in prestellar cores (Crapsi et al. 2007, Pagani et al. 2007) and dark globules (Pineda 
& Bensch 2007). The values for the cosmic-ray heating energy per ion pair available in the 

l ln this and other papers, the cosmic-ray ionization rate (h per H nucleus is used, which 
is smaller by a factor of two than the rate Ch 2 P er H 2 molecule of Indriolo et al. (2007) . 
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literature range over a factor of three. In this paper we will attempt to clarify this situation 
for both cosmic-ray and X-ray heating. 

The interaction of fast electrons slowing down in molecular gas was first analyzed more 
than three decades ago by Glassgold & Langer (1973, hereafter GL73), Cravens et al. (1975) 
and Cravens & Dalgarno (1978). Although GL73 was the first paper to consider cosmic ray 
heating in interstellar molecular regions, it suffered from the incompletely known electron 
cross sections of the early 70s, a crude energy-loss calculation, and cumbersome notation. 
Cravens et al. (1975) did a better job on cross sections and energy loss, but restricted 
themselves to low-energy electrons in order to evaluate the validity of the commonly used 
continuous slowing down approximation. This limitation was removed by Cravens & Dal- 
garno (1978). There are no glaring discrepancies between Cravens & Dalgarno (1978) and 
GL73. All of these papers considered pure H 2 regions, i.e., the roles of H + and He + ions 
in molecular gas were ignored. A more up to date and complete analysis was carried out 
by Dalgarno et al. (1999, hereafter DYL). DYL considered carefully all of the energy loss 
channels for electron energies from 30 eV to 1 keV in various mixtures of H, H2 and He. 
They showed how the energy expended to make an ion pair W is partitioned between elastic 
and non-elastic processes, but they did not fully treat the heating. 

Cosmic-rays and X-rays (or the equivalent, fast electrons) produce ions and excited 
molecules that can interact with the dominant neutral atomic or molecular gas. The products 
of these reactions carry away a significant amount of the available energy and heat the gas. 
We refer to this as chemical heating, and quantify it by a quantity Qchcm defined per ion 
pair. Chemical heating by cosmic-rays and X-rays occurs in a wide range of applications from 
diffuse and dense interstellar clouds to pre-stellar cores, protoplanetary disks, and planetary 
atmospheres. In this paper we use the results of DYL to inquire how the excited species and 
ions heat a gas mixture of H, H 2 , and He. We will show that ~ 1/2 of the energy of X-rays 
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and cosmic rays can go into heating in dense molecular regions and that chemical heating 
can be the most important part of the heating. We also note that chemical heating occurs 
when molecular regions are irradiated by far ultra-violet (FUV) photons. 

This paper is organized as follows. In Section [2] we outline the physical basis of our 
analysis and make quantitative estimates in Sections [3] and HI We discuss the most important 
results in Section [5j The paper ends with a short summary in Section |6j 

2. The Physical Basis for Cosmic-ray and X-ray Heating 

A basic tenet of earlier work on cosmic-ray heating is that proton cross sections can 
be represented accurately by the corresponding ones for electrons (which are generally more 
available), at least at sufficiently large energy. Cravens & Dalgarno (1978) discussed this 
explicitly on the basis of the good agreement between experiments on proton collisional 
ionization of H 2 (Hooper et al. 1961) and the Born approximation, as developed by Bates & 
Griffing (1953), all the way down to a proton energy of 0.06 MeV. This conclusion stems from 
the fact that high-energy ionization cross sections are a function of the incident particle's 
speed, so proton cross sections are the same as those for electrons with the same speed, i.e., 
when the electron energy E satisfies the relation E = (m c /m p )E p . 

An energetic cosmic-ray proton loses energy in cosmic gases in small steps that are 
characterized by the energy W to make an ion pair, known to be close to 37 eV from extensive 
experimental and theoretical studies. In their interactions with the hydrogen and helium 
of the gas, each cosmic-ray (proton) ionization generates a distribution of relatively slow 
electrons with average energy of 19 eV (Opal et al. 1971; GL73) and that leads to another 
"secondary" ionization two-thirds of the time. X-rays with keV energies develop a similar 
sequence once the initial photon is absorbed by a heavy atom by L- or K-shell electron 
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ejection. The energy of this primary electron is E — hvx — B — £, where B is the binding 
energy of the K or L shell of the heavy atom and £ is the mean energy of the resulting ion. 
It is significantly less than the initial X-ray energy, but much of the difference is restored in 
the form of several prompt supra-thermal Auger electrons with keV energies. There is also 
some X-ray fluorescence that, by the processes just described, generates more fast electrons. 
The net effect of X-ray ionization of heavy atoms, accompanied by Auger electrons and 
fluorescence, is that most of the initial X-ray photon energy is recovered in the form of the 
primary electron plus several other fast electrons. Thus, once the the processes attending 
the absorption of a keV X-ray have occurred, the main ionization effects of the X-ray are 
accomplished by the fast electrons that it produces. This is almost exactly the situation 
with cosmic-ray ionization, except that these electrons are replaced by the cosmic-ray proton 
which, as noted above, behaves like an electron of the same speed. This is the basis for our 
using the extensive results of DYL on the interaction of electrons with atomic and molecular 
gases to calculate the heating by cosmic-rays and X-rays. 

The spectra of cosmic-rays and X-rays can change from region to region due to proximity 
to a source or absorption within the cloud of interest. One example is the absorption of 
low-energy protons in thick clouds, proposed by Indriolo et al. (2009) and Padovani et 
al. (2009) to explain the reduction of the ionization rate inside dense molecular clouds. 
Another is the observed variation in the observed spectra of the X-rays emitted by young 
stellar objects. DYL give results for mono-energetic electrons with energies from 30 to 
1000 eV. Their results do not vary by more than 10% for E > 200 eV, and if the underlying 
spectrum does not extend below this energy for electrons or the equivalent energy of 0.4 Mev 
for protons, the variation of the DYL results at low electron energy can be ignored. Otherwise 
the results presented in this paper have to be integrated over the appropriate spectrum. 
Below E = 30 eV (or ~ 60keV protons), the lowest energy considered by DYL, one can use 
GL73, but the equivalence between electrons and protons has very likely broken down at this 
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point. In this paper, we focus on electrons with energies greater than 500 eV, corresponding 
to protons with energies greater than 1 MeV. 

We use DYL to analyze all of the outcomes of a fast electron as it slows down by 
interacting with a cosmic gas-mixture dominated by H, H2, and He, ignoring the much 
less abundant heavy elements. Unlike DYL, we make a comprehensive effort to evaluate 
how much of the energy expended in each process goes into heating - the primary goal of 
this paper. For this purpose, the most useful results in DYL are Table 4 (ion production), 
Table 5 (excitation), and Table 7 (heating efficiency), all for the case of H, He and H2, He 
gas mixtures. 

The energy partitions for a 1-keV electron in a He, H 2 mixture in the limit of zero 
electron fraction are listed in Table 1. The entries are independent of energy above 100-200 
eV. If the X-ray spectrum extends below 200 eV or the cosmic-ray spectrum below 400 keV, 
the full tables in DYL contain the information needed for lower energies. Table 1 has been 
organized into four parts: ion production; excitation of H 2 levels; excitation of He levels; and 
the heating from elastic scattering combined with the heating from rotational excitation. All 
energies are expressed in eV. The first column names the process; the second column the 
energy expended for that process W; the third column the number of such processes for an 
electron energy E, N = E/W, whe the electron energy is E — 1 keV; the fourth column the 
threshold energy E th for the process; and lastly the energy associated with the process, i.e., 
the product NE th of the two previous columns (for E = 1 keV). 

For example the first four lines of Table 1 show that a keV electron generates 23.9 H^ 
ions, 1.1 H + ions, 2.2 He + ions, and < 0.1 He 2+ ions, for a total of 27.2 ions. This follows 
from the fact that an incident electron expends on average an energy Wi to make an ion 
pair of type i (i — 1 for Hj, % = 2 for H + , % — 3 for He + , and % = 4 for He 2+ ). More 
generally, if Ni = E/Wi is the number of ions generated of type i, the total number of ions 
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Table 1. DYL Energy Partitions for a He, H 2 mixture 



Process 


W 


N 




NE th 




(eV) 




(cV) 


(eV) 




41.9 


23.87 


15.44 


369 


H+ 


921 


1.09 


18.08 


20 


Hc+ 


459 


2.18 


24.59 


54 


Hc 2 + 


16000 


0.0625 


54.42 


3 


Total Ion production 








446 


B 


117 


8.55 


11.37 


97 


C 


132 


7.58 


12.42 


94 


Dissociation 15 


92.6 


10.8 


~ 13 


140 


H+H(2p) 


534 


1.87 


14.7 


28 


v = 1 


7.81 


128 


0.516 


66 


v = 2 


109 


9.17 


1.032 


10 


Total H2 excitation 








435 


He 2*S c 


15000 


0.0667 


20.62 


1 


2 X P 


1940 


0.515 


21.22 


11 


2 3 S 


34900 


0.0287 


19.82 


1 


2 3 P + n > 2 3 


22500 


0.0444 


21.0 


1 


71 > 2 1 


3520 


0.284 


23 


7 


Total He excitation 








21 


Elastic collisions and 










Rotational cxcitation d 








57 


Grand total 








959 



a For the case x e = and E = 1 keV 
b Dissociation occurs mainly via triplet excitations. 
c He excitation follows the notation in DYL. 
d Based on DYL "heating efficiency" in their Table 7. 
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is N — J2 i Ni = E J2i W { 1 and the average energy per ion pair is 

E 



CC^- 1 )- 1 - 



N 

According to Table 1, E/N = 36.8 eV per ion pair for E — 1 keV and zero electron fraction. 

The data in DYL account for 96% of the initial electron energy E — 1 keV in Table 
1. The small deviation from 100% is not significant in light of the uncertainties in the 
calculations. It may be due to the use of the approximate fitting formulae for the parameters 
W, for which DYL quote an accuracy of 5-15%. Another possibility is the omission of triplet 
states starting near 14 eV that lie above the a, b, and c levels, which DYL refer to as missing 
"pseudo states". In addition, most of the cross sections used in the energy loss process are 
probably not known to an accuracy of 10%. Approximately 47% of the energy in Table 1 is 
accounted for by ionization, 47% by excitation, and 6% by direct heating and the collisional 
de-excitation of H2 rotational levels. Applying this to the 36.8 eV to make an ion pair, elastic 
scattering and rotational excitation lead to only about 2 eV heating per ion pair in neutral 



gas. However, as the ionization fraction increases beyond x e 



, the direct heating 



increases, largely due to Coulomb collisions with ambient electrons 2 ]. Most cosmic-ray and 
X-ray heating in molecular regions comes from the reactive ions and excited states. This is 
exactly the part of the energy loss problem that is not treated by DYL. 



3. The Heating Energy 

In general the heating energy per ion pair Q consists of contributions from collisions 
with H and H 2 , 

_ x(H)Q(H)+x(H 2 )Q(H 2 ) 
* x(H) + x(H 2 ) ' [A) 



2 It is important to note that the DYL considerations do not apply above x e = 0.1. 
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where x(H) + 2£(n 2 ) = 1; the functions Q(H) and Q(H 2 ) depend on the abundance of He 
and especially on the abundance of electrons. The volumetric rate of ion-pair production is 
Ch^H; where Ch is the ionization rate per H nucleus and tt-h is the volumetric density of H 
nuclei. Thus the volumetric heating rate is given by T = Q Ch^h- 

The dominant heating process in atomic regions is elastic scattering. DYL include 
Coulomb collisions with ambient electrons, but only for the case where the electron abun- 
dance relative to total H nuclei is < 0.1. For negligible ionization, it is twice that for H2 
regions, as can be seen in Table A of the Appendix, the analog of Table 1 for atomic regions. 
It is included in Q(H), the first term in the numerator of Eq. |5J In Appendix A we show 
that chemical heating induced by the atomic ions H + and He + reacting with H is very small. 
As shown in DYL Figure 8, the dominant excitation product in atomic H is the 2p level of 
atomic hydrogen. The resultant Ly-a photons will either escape or be absorbed by and heat 
the dust, and not the gas. 

The heating in molecular regions is the sum of the effects from elastic collisions plus 
rotational excitation (Q e i/rot); excitation of H 2 vibrational levels (Q v ib) ; dissociation of H 2 
(Qdiss), and chemical heating (Q c h em ) : 

Q(H 2 ) = Qcl/rot + Qvib + Qdiss + Qchem ■ (3) 

In the following paragraphs we discuss the first three contributions to <5(H 2 ), and then focus 
on Qchem in the next section. 

3.1. Elastic Scattering and Rotational Excitation 

The critical densities for the rotational levels of H 2 depend on the transition and also 
on the H/H 2 ratio and the temperature (e.g., Le Bourlot et al. 1999). In many situations the 
lowest S(0) and S(l) rotational levels are collisionally de-excited and most of the rotational 
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excitation goes into heating. GL73 estimated that a density n-^ = 1000 cm -3 would satisfy 
this condition on density, independent of temperature. We use DYL for the combined heating 
due to elastic scattering and rotational excitation, 

Qei/rot(H 2 , e) = Q d (H 2 , e) + Q rot (H 2 , e) . (4) 

This quantity is then expressed in terms of a "heating efficiency" 77 for a H2, He gas mixture, 

g e i/rot(H 2 ,e) = ^(H 2 ,e)W(H 2 ,e), (5) 

where W r (H 2 , e) is the average energy to make an ion pair. Both £h and the heating function 
Q are defined in terms of ion pairs. For the H, He gas, there is of course no rotational 
heating, and the heating is expressed in a similar form where the heating efficiency ?y(H, e) 
now describes only elastic collisions, 

Q ol (H,e)=7 ? (H,e)W(H,e). (6) 

However, W(H,e) ~ W(B. 2 ,e) to an accuracy of ~ 5% for electron energies greater than 
several hundred eV, and the first term of Eq. [3] becomes, 

n m \>n m ^ ~ [x(H)y(H,e) + x(H 2 )^(H 2 , e)] 

Qoi(H, e) + <2ei/rot(H 2 , e) w x(H) +x(H 2 ) W ( R ^ e ) ( 7 ) 

with W(H 2 ,e) w 37 eV asymptotically for E > 500 eV. 

The heating efficiencies r] are expressed by the fitting formulae, DYL Eq. (14), 

V= I - — — . (8) 

The two sets of fit parameters are given at the bottom of DYL Table 7 ("Two-Gas Mixtures" ) 
as a function of electron energy E. For example, the E — 1 keV values are, 

7?(H,e)o = 0.117, a = 0.678, C = 7.95, (9) 
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for atomic regions, and 

r/(H 2 ,e) = 0.055, a = 0.366, C = 2.17, (10) 

for molecular regions. These values might be appropriate for the keV X-rays emitted by 
young stellar objects. These fits have only a modest accuracy ~ 10 — 15%. 

According to DYL (see Table 1), Qei/rot accounts for only 57 eV, or ~ 6% of the energy 
of the 1-keV incident electron. The rest of the energy is almost equally shared by ionization 
and excitation of He and H 2 . Thus, the DYL heating efficiency for molecular regions is only 
a small part of the total. With W ~ 37 eV to make an ion pair, elastic scattering and 
rotational excitation lead to Q e \/mt = 2.1 eV heating per ion pair in neutral H 2 , He gas. This 
is half of the direct heating in an atomic H, He, ~ 4.3 eV under the same conditions. 



3.2. Dissociation Heating 

The most important pathway to dissociation is collisional excitation of H 2 triplet states 
that start at 11.9 eV (a and c), 13.4eV (e) and possibly above. These levels decay into 
continuum states b of the triplet repulsive potential and have essentially unit probability 
for dissociation. Unlike the Lyman and Werner transitions to the singlet B and C levels, 
these spin-flip transitions are forbidden for photon excitation from the singlet ground state 
X. Table 1 indicates that for E — 1 keV and x e = 0, there are 10.8 dissociations per keV. 
Adopting the DYL estimate of 5.4 eV per dissociation, we obtain 58 eV heating per keV, or 
2.14eV per ion pair for a neutral H 2 , He mixture []. More generally we use DYL Table 5 for 

3 This number does not include the heating from dissociation following collisional- 
excitation of the B and C states, which is much smaller, only ~ 0.6 eV per keV. 
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the dissociation heating in Eq. [3J, 



x(E 2 ) D 



(11) 



diss 



x(H) + z(H 2 ) 1 + Cx% 



where D { 



'o = 



2.14 eV, a = 0.574 and C = 22.0. 



3.3. Vibrational Heating 



Fast electrons can excite the vibrational levels of H 2 by direct collisions or by fluorescent 
de-excitation of electronically excited B and C levels. The excitation energy goes into 
heating if the densities are high enough for the levels to be collisionally de-excited. Recent 
calculations of the collisional de-excitation rate coefficients (e.g. Wrathmall et al. 2007 for 
H; Lee et al. 2008 for H2; Balakrishnan et al. 1999 for He) replace earlier work of Tine 
et al. (1997) and Le Bourlot et al. (1999). The coefficients are generally much larger for 
collisions with atomic H than for H2 and He. They are very small below 100 K and increase 
rapidly above this temperature. Thus one can expect little quenching of vibrational levels in 
the cool and moderately dense regions of molecular clouds. By contrast, in the surface layers 
of the inner regions of protoplanetary disks, where nn ~ 10 10 cm~ 3 and T ~ 1000 K, the 
densities of both H and H 2 are large enough to collisionally quench vibrationally-excited H 2 
molecules and provide significant heating. In less dense molecular regions, atomic H collisions 
may be effective in collisionally de-exciting H 2 , but this depends on the temperature and the 
H/H 2 ratio as well as the total density %. 

In order to obtain the heating from direct collisional excitation in dense H 2 regions, we 
use DYL Table 5 for H 2 , He mixtures which gives the following fit for the energy to excite 
the v = 1, 2 levels, 



The parameters for E = 1 keV are given in Table 2 for E > 500 eV, with W v and AE V (the 



e v = W v (l + C v x?). 



(12) 
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vibrational excitation energies) in eV. 

Because 1000/et, is the number of excitations per keV of electron energy, the high-density 
direct vibrational heating per keV of incident electron energy is, 

1000 . *\ „ 

AE X + AE 2 eV . 



/1000 . „ 1000 . „ \ 



The vibrational heating per ion pair is obtained by dividing this expression by the number 
of ion pairs per keV (27.2 for E > 500 eV), with the result, 



Qdir/vib^ * ( " 2) 19.0eV 



1 2 
— + — 

ei e 2 



(13) 



In a pure H 2 region with negligible ionization, the vibrational heating per ion pair is 2.8 eV 
for E > 500 eV, with v = 2 contributing about 15%. 

Excited vibrational levels are also produced when collisionally excited B and C states 
decay with the emission of fluorescent photons in the 1500-1600 A range, with a maxi- 
mum near 1575 A. We can calculate the high-density heating that follows the collisional 
de-excitation of the vibrational levels using the above method for direct collisional excita- 
tion based on Eq. [T2J The necessary parameters again come from DYL Table 5 for H 2 , He 
mixtures. 

Table 2. Direct Vibrational Parameters 8, 



v W v (eV) C v a v AE V (eV) 



1 7.81 23,500 0.955 0.516 

2 109 10,700 0.907 1.032 



E = 1 keV 
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With roughly 4eV in vibrational excitation energy and 27.2 ion pairs per keV, the 
high-density B, C vibrational heating per keV of incident electron energy is, 



1000/ tB + lOOOte 
27.2 



1 1 



The B, C vibrational heating per ion pair is then 

z(H 2 ) 



QBC/vib — 



147 eV 



1 1 

es e c 



(14) 



x(H) + x(H 2 ) 

If we substitute the values in Table 3 for x e = 0, the result is 2.4 eV per ion pair for a neutral 
H 2 , He gas, a small but not negligible contribution to the high-density heating. 



As discussed earlier, both Eqs. [13] and [14] hold only if the density well exceeds the 
critical density n CT for the de-excitation of the vibrational transitions of the H 2 ground level. 
They therefore should be multiplied by an appropriate factor O vi b, that depends on the ratio 
n cr/ n H, such that vi b = for small n cr /nH and v ;b = 1 for large n cr /nn. Their contribution 
to the total heating in Eq. [3] is now, 



Qvib — ©vib (Qdir/vib + QBC/vib) ■ 



(15) 



The B and C levels have a ~ 15% probability to dissociate (e.g., Abgrall et al. 1997), 
with a yield of approximately 0.25 eV in typical interstellar conditions (Tielens 2000). The 



Table 3. B, C Vibrational Parameters 3 



Level W v (eV) C v a 



B 117 7.09 0.779 

C 132 6.88 0.802 



E = 1 keV 
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small branching ratio for dissociation and the small energy yield mean that the heating per 
ion pair following dissociation is very small, e.g., 0.02 eV for E = 1 keV, and thus negligible. 

4. Chemical Heating 

Chemical heating derives mainly from reactions instigated by the primary cosmic-ray 
and X-ray ions, H^, H + , and He + , with neutral species and electrons. It was first considered 
in the context of EUV heating of Jupiter's upper atmosphere (Henry & McElroy 1969) and 
then for cosmic ray ionization of interstellar molecular clouds (GL73). These authors focused 
on H^~, its transformation into Hjj", and the destruction of by dissociative recombination. 
In dense molecular regions, however, where x e <C ^(H 2 ), dissociative recombination has 
to compete against ionic reactions with neutral species. Rather than treat all of the ion- 
neutral reactions that can be traced to the primary ions, H^", H + , and He + , we focus on only 
the potentially most abundant neutrals: CO, H 2 0, and O. Although this treatment is not 
completely general, it should suffice to demonstrate the nature and magnitude of chemical 
heating in molecular regions. The chemical heating will be expressed in terms of its values 
for the three main initiating ions in H2 regions, 

Qchem(H 2 ) = Qchem (H 2 , H 2 ) + Q c hem(H + , H 2 ) + Q c h em (He + , H 2 ) . (16) 

According to the discussion of Table 1 in Section El each of these terms includes a weighting 
factor equal to the fraction for each ion of the total, i^Hj) = 0.88, F(H + ) = 0.04, and 
F(He + ) = 0.08, ignoring the tiny fraction of He 2+ . In Appendix A we show that chemical 
heating in atomic regions is negligible. 
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4.1. Chemical Heating of H 2 

As discussed in Section 2, the most abundant ion generated by cosmic-rays and X-rays 
in H 2 regions is . We first consider its main destruction routes; H + and He + ions will be 
analyzed later in this section. is mainly destroyed by dissociative recombination and by 
proton transfer with H 2 , 

H+ + e^H + H, H+ + H 2 ^H+ + H. (17) 

with rate coefficients, (3 = 2.0 x 1CT 7 T _1 / 2 cm 3 s _1 (a rough fit to Schneider et al. 1994) and 
k = 2 x 10" 9 cm 3 s _1 (Theard & Huntress 1974). The probability for H3 production is, 

pm + tj+a _ MH 2 ) x(H 2 ) 
1 2 ' 3j kx(R 2 ) + f3x c x(H 2 ) + lOOT' 1 ^ ' ^ 

so that if, 

x e < 0.1 |x(H 2 ) = 0.01 (— ^) x(H 2 ) , (19) 

the rate of dissociative recombination is less than 10% that for H 2 destruction. For regions 
with considerable H 2 , this condition is well satisfied, P^^fht) 1, and most of the H 2 
ions are transformed into Hg\ 

The destruction pathways for H3 are first, dissociative recombination with two branches, 

H+ + e ->• H 2 + H (25%) and H + H + H (75%) , (20) 

and total rate coefficient (3' = 4.5 x 10~ 6 T _a65 cm 3 s _1 (Sundstrom et al. 1994), and second, 
proton-transfer reactions such as, 

H+ + CO ->• HCO+ + H 2 h = 1.6 x 10~ 9 cm^" 1 
H+ + H 2 ->■ H 3 0+ + H fc 2 = 5.3 x 10~ 9 cmV 1 (21) 
H+ + O ->■ OH+ + H 2 k 3 = 0.8 x 10~ 9 cmV 1 . 
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Dissociative recombination, Eq. [201 was considered by GL73. It is competitive with the 
reactions in Eq. |5T] where the electron fraction is relatively large. Denoting the abundance 
of the neutral reactants in Eq. [2D by X{ with i = 1, 2, 3 for CO, H 2 0, and O, respectively, we 
introduce branching ratios for the main reaction pathways of H^, 

Bc= yk P ' X ° t B *= v h klX lR> (* = 1 > 2 ' 3 )- ( 22 ) 

Zjj/CjXj -+- p X e Z-nKiXi -\- p X e 

In order for dissociative recombination to compete with ionic reactions at the 10% level, the 
electron fraction must satisfy the condition, 

x e > 0.1 = 2.22 x 10- 9 T ' 65 (EiksO-ia , (23) 

where (£jfcja;j)_i3 = (Sj/cjXj/10~ 13 cm 3 s _1 ) measures the reactivity of hht with abundant 
neutrals in units determined by neutral abundances of order 10 -4 and ionic reactions with 
rate coefficients of order 10~ 9 cm 3 s _1 . Even in cold molecular regions, (T ~ 10 K) dissociative 
recombination can still play a role in the presence of neutral reactions. Eq. [23] for H^" is much 
less restrictive than Eq. [19] for H^ because of the reduced abundance of heavy elements 
compared to hydrogen. 

The first reaction in Eq. [21] with CO is simple to treat because the product HCO + 
is mainly destroyed in one step by dissociative recombination back to CO + H, whereas 
dissociative recombination of H 3 + in the second reaction has three branches, H 2 + H, 
OH + H2, and OH + 2H. The situation for the reaction with O is similar because OH + in 
the third equation is quickly transformed into H3O 4 " by hydrogenation reactions with H 2 . 
In this sequence, we can ignore dissociative recombination of OH + and H 2 + because of 
conditions similar to Eq. [19] 

Keeping in mind that the probability for the conversion of H^ into Hjj", P(H 2 f ,H^), is 
essentially unity according to Eqs. [18] and [191 the chemical energy released by the creation 
of an H^ ion can be obtained by consolidating the reactions in Eq. [17] and Eq. [21] into the 
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equivalent reactions, 

e + H+ -> H + H (24) 

e + H+ + H 2 + H 2 -> H 2 + H 2 + 2H, OH + 2H 2 + H, OH + H 2 + 3H (25) 

e + H+ + + 2H 2 -> H 2 + 4H, OH + H 2 + 3H , OH + 5H (26) 

where q x = 11.1 eV, q 2 = 7.8 eV, and q 3 = 5.7eV are the net energy yields for each of 
the three consolidated reactions of Hj, Eqs. [MJ [251 and [2S The g are averages over the 
outcomes generated by the three branches in the dissociative recombination of H 3 + , 

H 3 0+ + e -> H 2 + H (25%); OH + H 2 (15%); OH + 2H (60%) . (27) 

For example, Eq. [21] results from adding the full sequence of reactions that are involved when 
H^ reacts with CO: 

H+ + H 2 -> H+ + H 
H+ + CO -> HCO+ + H 2 
e + H+ -> CO + H. 

The chemical heating due to the reactions of the H^ ion is then the rate at which that 
ion is produced per unit volume multiplied by the probability that it is transformed into an 
H^ ion, (Eq. I18p . times the heating averaged over the three branches of the latter's reactions 
and including the heating from the dissociative recombination of H^, all multiplied by the 
H 2 fraction. In terms of the notation of Eq. [161 the chemical heating per ion pair stemming 
from the production of the H^ ion is, 

Qc hcm (H+) = -^^Lf(H 2 + )P(H+H 3 + ) p i S j (H+)g i (H+) +B e g 6 (H+)] . (28) 

The heating energies q { are given following Eq. [261 an d the branchings are defined in Eq. [221 
g c (H 2 f ) = 7.6 eV. For the case where ion-molecule destruction prevails over dissociative re- 
combination, we can estimate the value of the chemical heating for the case where the 
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abundances for CO, H 2 0, and O all equal 10 4 . The result is Ej 5j(H^") g«(H^) = 8.4 eV, 
but this must still be multiplied by F(H+) = 0.88 to yield Q ch em(H^) = 7.2 eV. 

Examination of Eqs. 1251 and [261 reveals that the products of the ion-molecule reactions 
for H2O and O involve the radicals OH and H whose further reaction can lead to more 
chemical heating. For example, OH can be converted to H 2 by the exothermic reaction, 

OH + H 2 -)• H 2 + H, (29) 

and the H atoms can make H 2 by formation on dust grain surfaces "gr" , 

H + H + gr H 2 + gr'. (30) 

The energy yield of the OH reaction is 0.65 eV. When this is added to the energy yields of 
Eqs. [25] and [261 the chemical heating produced in these reactions is increased slightly, q 2 = 
8.4 eV and q 3 = 6.4 eV. The previous estimate of T li B i q i for the case of equal abundances 
for CO, H 2 0, and O (all 10" 4 ) is increased from 8.4 eV to 8.8 eV, and Q chein (H+) = 7.7 eV. 

The heating from Eq. [30] depends on how much if any of the energy release goes into 
kinetic energy of the newly formed molecule compared to internal excitation and excitation 
of the birth grain. Theoretical estimates and guesses in the literature range from 0.5 to 
3.0 eV, but laboratory experiments (e.g., Roser et al. 2003) suggest that the H 2 molecule 
thermally accommodates to the temperature of the grain before making its final escape, at 
least as far as kinetic energy is concerned. In a recent experiment, Lemaire et al. (2010) 
present evidence that ~ 30% of newly formed H 2 molecules are vibrationally excited for 
dust temperatures as warm as 70 K. Assuming that the results of Lemaire et al. applies to 
molecular regions where the density is large enough for vibrationally excited H 2 molecules 
to be collisionally de-excited, the chemical heating for the H 2 and O channels would be 
increased by 1.5 eV. This would then lead to an increase in the heating function by about 
leV to Qchem(H 2 f ) ~ 9 eV. In the numerical estimates below, we ignore this contribution 
pending the resolution of this long standing issue regarding H 2 formation on grains. 
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4.2. Chemical Heating of H + and He + 

In principle, the chemical heating of H + and He + ions can be treated in a similar way 
as H^, but the outcomes depend on the chemical composition of the molecular region in 
question. Thus H + does not react with CO, the most abundant and stable heavy molecule, 
and its fast charge-exchange with O is closely balanced by the reverse reaction. It does react 
with H2O, which can be very abundant under certain circumstances, and with many organic 
species that usually have small abundances. In a cold molecular cloud, where H 2 may be 
expected to be frozen out on grains, this route will be shut down. On the other hand, in 
a dense warm region like those observed in the inner regions of protoplanetary disks, the 
primary destruction pathway for H + is 

H+ + H 2 -)• H 2 0+ + H, (31) 

followed by hydrogenation to H 3 + , 

H 2 + + H 2 -)• H 3 + + H, (32) 

and destruction of H30 + by dissociative recombination, according to Eq. [27] The energy 
balance equation is similar to Eq. [25] for reacting with H 2 0. The only difference is that 
is replaced by H + on the left side and H 2 is replaced by H on the right side. The net 
change in exothermic reaction yield (and in the maximum heating) is just the difference 
between the ionization potentials of H 2 and H, or 1.83 eV. Thus the chemical heating arising 
from H + analogous to Eq. [28] is, 

a hcm (H + ) = 2 ) F(H+) B 2 (H+) g t (H + ) , (33) 

where F(H + ) is the fraction of H + ions (0.040 from Table 1), 5 2 (H + ) ~ 1 is the branching 
ratio for H + to react with H 2 0, and g 2 (H + ) = ^(H^ ) — 1.83 eV. In principle, B 2 (R + ) is given 
by an expression similar to Eq. [22] with 0', the dissociative recombination rate coefficient 
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for H3 , replaced by the radiative recombination rate for H + . Because the latter is smaller 
by many orders of magnitude, destruction of H + by radiative recombination can be ignored. 
For a completely molecular region, F(H. + ) = 0.04, B 2 ~ 1, and g 2 (H + ) = 6.9 eV, resulting in 
Qchem(H + ) = 0.28 eV. This estimate includes the energy that can be recovered from the OH 
radical but not from the conversion of product H atoms to form H 2 on grains. It is much 
smaller than the heating per ion pair for H^" because H + is less energetic and because only 
4% of the ion pairs involve H + . This estimate is also sensitive to the physical conditions 
because the H 2 abundance can vary, especially due to freeze out on grains, and because 
reactions with other species may occur such as charge transfer with neutral heavy atoms. 

The heating by the He + ions is more complicated because it reacts with many molecules 
including H 2 . Here we consider the dominant reaction to be, 

He+ + CO C+ + + He, (34) 

followed by, 

C + + H 2 -> HCO+ + H, (35) 

and then by dissociative recombination of HCO + . A similar analysis to that followed above 
for H2 gives an energy budget equation, 

e + He + + H 2 -> O + 2H + He , (36) 

with an energy yield of 15.0 eV. Including the reactions of the remaining O atom gives 
another 0.6 eV, or gi(He + ) = 15.6 eV. On this simplified basis, the heating due to He + ions 
is 

Q chcm (He + ) = x(H ^ (H2) ^(He + ) Bi(He+) 9l (He + ) . (37) 

Using F(He + ) = 0.80, i?i(He + ) = 1 (radiative recombination is again unimportant in molec- 
ular regions), Q c hem(He + ) = 1.25 eV. This may be a slight overestimate because other reac- 
tions of He + have been neglected, e.g., with water. 
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4.3. Summary of Chemical Heating 

To summarize this treatment of chemical heating, we have considered the energy that is 
available for heating due to the reactions of , H + , and He + ions, which account respectively 
for 88%, 4%, and 8% of the ~ 27 ion pairs per keV produced by high energy electrons 
(E > 500 eV). The ion accounts for 7.7eV per ion pair, which is most of the chemical 
heating (84%). This is very close to the heating from the dissociative recombination of hLt , 
~ 7.6 eV (GL73). Even the total chemical heating, 9.2 eV per ion pair, is insensitive to 
the electron fraction because the chemical heating from H + and He + are also essentially 
independent of the electron fraction. There are a number of fine points that could add or 
subtract 1-2 eV from this result, depending on the physical conditions and on the role of 
poorly understood processes. For example, this treatment of high-density chemical heating is 
an over-estimate in that no energy is assumed to be lost by vibrationally-excited molecules 
produced in the relevant chemical reactions, on the assumption that the density is high 
enough for collisional de-excitation to be effective. Similarly, we have ignored the uncertain 
possibility that ~ 1.5 eV of excitation energy of newly formed H 2 molecules may be available 
for heating at high densities. 

5. Results 

In the previous sections we developed a theory for the heating of molecular gases exposed 
to cosmic rays and X-rays, with detailed consideration of chemical heating. To illustrate the 
magnitude of the various contributions to the heating, we list in Table 4 the maximum 
values for the case of zero electron fraction and densities high enough for vibrationally- 
excited molecules to be collisionally de-excited; the electron energy is E — 1 keV. The 
chemical heating in Table 4 is based on somewhat arbitrary choices for the most abundant 
neutrals: x(CO) = x(0) = rr(H 2 0) = 10~ 4 . Starting at the bottom, we see that the total 
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maximum heating is Q = 18.7 eV per ion pair; close to 50% of the energy expended per ion 
pair, which in this case isW = 36.8 eV. The largest fraction, about one half, is accounted for 
by chemical heating (9.3 eV). This value is ~ 20% larger than the value of 7.6 eV obtained 
from reactions Eq. [TH H£ + H 2 -> + H, and Eq. [2UJ H.[ + e — > H 2 + H, along the lines of 
GL73. The next largest contribution, 5.2 eV, comes from vibrationally excited H2 molecules, 
produced by direct collisional excitation or following the excitation of the B and C states. 
Elastic scattering, rotational excitation, and dissociation account for the rest, or ~ 20%. 

The values in Table 4 are based on several assumptions: the electron fraction is very 
small, the density very high and the abundances of C, H 2 and O all equal to 10 -4 . The 
chemical heating appears to be insensitive to the electron fraction as long as Eq. [T9] is 
satisfied. On the other hand, Coulomb collisions with ambient electrons begin to play a 
significant role in heating wherever x e approaches the 10 -4 level. This effect is automatically 
included in Eq. [71 The entries for vibrational heating in Table 4 are the most sensitive to 
density because very high densities are required to collisionally quench vibrationally excited 
H 2 molecules, as discussed in Section 2.4. This restriction means that vibrational heating 
is ineffective in molecular clouds outside of pre-stellar cores. The chemical heating also 
depends weakly on the abundances of the neutral species. For example, if instead of the 
abundances used for the estimates in Table 4, x(CO) = x(0) = x(H 2 0) = 10~ 4 , we choose 
x(CO) = 10~ 4 , x(0) = 2 x 10~ 4 , and x(H 2 0) = 0, the chemical heating from increases by 
about ~ 0.2 eV. Thus the chemical heating associated with the ion is reasonably robust 
in molecular regions at ~ 8 — 9 eV. 

Under certain circumstances, such as very low energy cosmic-rays or soft X-rays, the 
considerations of this paper need to be extended to electron energies E < 1 keV. Without 
attempting to treat this subject in detail, we can obtain a preliminary view of the situation 
from Table 5, which gives heating functions Q and the energy to make an ion pair as a 
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Table 4. Heating per Ion Pair a 



Process 


Q (eV) 


Q el/rot 


2.1 


Qdiss 


2.1 


Qdir/vib 


2.8 


<?BC/vib 


2.4 


Qchem(H^", H 2 ) 


7.7 


Qchem(H + , H 2 ) 


0.3 


Q chcm (H+,He + ) 


1.3 


Total heating Q 


18.7 



a E = 1 keV 
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function of electron energy. Many of the heating functions increase with decreasing energy, 
especially collisional excitation and dissociation. However the percentage of the energy to 
make an ion pair that goes into heating remains close to 50% down to E — 50 eV. Indeed, 
the following quantities vary by less than 10% for energies greater than 200 eV: total heating 
Q, energy to make an ion pair W and the percentage heating. 

We can make the implications of the theory in this paper more concrete by estimating 
the heating for representative examples of interstellar and circumstellar matter that have 
significant amounts of H 2 . The estimates in Table 6 are based on the assumption that the 
electron energy is greater than 0.5 keV, corresponding to cosmic-ray proton energies greater 
than 1 MeV. For the well-observed line of sight towards ( Per (Ay = 1) (an early focus of the 
Hjj" observations, Indriolo et al. 2007), we use the modeling results for many observed species 
by Shaw et al. (2008), who fit the cosmic-ray ionization rate per H 2 with ( 2 ~8x 10 _16 s _1 . 
The densities in this diffuse cloud are too low for vibrational heating. The estimate for 
elastic/rotational heating may be too high because the density and temperature in this 
cloud may not be large enough for significant quenching of rotational excitation, according 
to the discussion in Sec. 3.1. The chemical heating in this case comes from dissociative 
recombination because, according to Eq. [231 it dominates over ion reactions. All of the 
heating estimates for (Per have been reduced by 70% to take into account that atomic H 
has not been fully converted into H 2 . The entries under Hjj~ destruction in Table 6 indicate 
whether the chemical heating is dominated by dissociative recombination of ( "DR" ) or 
ionic reactions ("I"). 

For the case of a molecular cloud clump, dissociative recombination plays a role in the 
destruction of Hjj~, but ion- molecule reactions are probably more important. Again, clumps 
are not dense enough for vibrational heating, and they also may not be warm enough to 
quench rotational excitation. For the inner region of a pre-stellar core, the density may 
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Table 5. Heating Q (in eV) per Ion Pair vs. Electron Energy 



Electron Energy E (eV) 


30 


50 


100 


200 


500 


1000 


Qe\/roi 


7.7 


4.4 


3.0 


2.5 


2.1 


2.1 


Qdiss 


12.0 


5.8 


3.9 


2.7 


2.3 


2.1 


Qdir/vib 


10.9 


6.1 


4.0 


3.3 


2.9 


2.8 


QBC/vib 


4.3 


3.3 


2.9 


2.6 


2.4 


2.4 


Qchem(H^", H2) 


8.6 


8.3 


8.0 


7.8 


7.7 


7.7 


<5chem(H + , H 2 ) 





0.1 


0.3 


0.3 


0.3 


0.3 


Qchem(H^, He ) 


0.2 


0.5 


0.8 


1.0 


1.2 


1.3 


Total heating Q 


43.7 


28.5 


22.9 


20.2 


18.9 


18.7 


Energy to make an ion pair W 


76.8 


52.9 


43.1 


39.6 


37.4 


36.8 


Number of ion pairs N = E/W 


0.39 


0.95 


2.3 


5.1 


13.4 


27.2 


Percentage heating 


57 


54 


53 


51 


51 


51 
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be high enough for vibrational heating, but their low temperatures suggest that the elas- 
tic/rotational heating may be somewhat smaller than given in Table 6. The freeze-out of 
volatiles, including some CO, means that the chemical heating will probably be close to the 
value for dissociative recombination. 

The densities in the inner molecular layer of a protoplanetary disk are often in excess 
of 10 9 cm -3 , and both rotational and vibrational heating should be effective. Because these 
regions are close to a stellar X-ray source, the X-ray ionization rates are much larger than 
the cosmic-ray rates for interstellar matter. Thus the electron fractions are relatively high 
at the top of the molecular region, ~ 10~ 7 — 10~ 6 , and Eq. [23] indicates that dissociative 
recombination still plays a role in the chemical heating of the very dense inner regions. 
However, for sufficiently large vertical columns, the X-ray ionization is sufficiently reduced 
that ionic reactions dominate chemical heating. 

All of the cases considered in Table 6 contain a significant if not an overwhelming 
fraction of H2. These rough estimates for the total heating per ion pair increase with the 
total density, and range from ~ 10 eV per ion pair for a moderately thick diffuse cloud like 
(Per to ~ 18 eV per ion pair for the densest regions. Perhaps the most important conclusion 
from these four examples is that cosmic-ray and X-ray heating are sensitive to the physical 
conditions, as expressed in the equations presented in Sections 3 and 4. The values of Q 
listed in Table 6 for various environments, combined with prescriptions for £h as function 
of the column density of the ambient gas (Padovani et al. 2009) and for the magnetic field 
intensity (Padovani & Galli 2011), allow the determination of the cosmic-ray heating rate in 
clouds, cores, and protostellar disks. 
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Table 6. X-ray and Cosmic-ray Heating (in eV) in Molecular Regions 





(Per 


molecular cloud 


prestellar core 


protoplanetary disk 




diffuse cloud 


clump 


inner region 


active region at 1 AU 


n H (cm 3 ) 


80 


10 4 


10 7 


10 10 


T(K) 


~ 60 


10 


6 


1000 


x e 


2 x 10~ 4 


10- 7 


10" 9 


lO" 6 


H3 destruction 


DR* 


dr + r 


DR + I 


DR + I 


<5ci/rot (eV) 


4 


2 


2 


2 


Qvib (eV) 








5 


5 


Qdiss (eV) 


1 


2 


2 


2 


Qchem (eV) 


5 


9 


8 


9 


Total heating Q (eV) 


10 


13 


17 


18 



*DR stands for dissociative recombination and I for ionic reactions. 



-30 - 
6. Conclusion 

Dalgarno et al. (1999; DYL) made an extensive study of the energy loss of fast electrons 
in H, He and H 2 , He gas mixtures. The electrons are at the heart of the interaction of 
cosmic rays and X-rays with interstellar and circumstellar matter. Although DYL analyzed 
essentially all excitation and ionization processes, they only discussed the heating from elastic 
collisions and rotational excitation of H 2 . Starting from the results in DYL, we have extended 
the scope of cosmic-ray and X-ray heating to include all of the relevant interactions. One of 
the main conclusions of this study is that heating by fast-electrons depends on the physical 
properties of the gas, i.e., on the abundance of electrons, H 2 molecules, and heavy atoms and 
molecules, and also on the total density of hydrogen nuclei. The electron fraction is important 
because, once it exceeds a certain level, heating by collisions with ambient electrons becomes 
important. The electron abundance also determines whether the destruction of the Hjj~ 
ions proceeds by dissociative recombination or ionic reactions, which affects the quantitative 
amount of chemical heating. Of course the H 2 /H abundance ratio is important because the 
diversity of the energy levels of H 2 offers many more channels for energy loss than atomic 
H. The dependence on physical conditions means that X-ray and cosmic-ray heating cannot 
be specified by a single number. This should be clear from the fact that in neutral atomic 
regions the heating efficiency is only 12%, whereas in neutral molecular regions it can reach 
50% at very high densities. 

A wide range of values for the heating per ion pair can be found in the literature on 
molecular clouds, from 7 eV (Stahler & Palla 2004), to 20 eV (Goldsmith & Langer 1978, 
Goldsmith 2001), but not necessarily for the reasons given here, e.g., in the above discussion 
of Table 6. Many authors adopt Goldsmith's value while others estimate intermediate values 
of 10-15 eV based on GL73, e.g., Maloney et al. (1996), Yusef-Zadeh et al. (2007), and 
Krumholz et al. (2011). As shown in Table 4, 13 eV is a good choice for not too dense 
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molecular gas. 

Chemical heating also applies to regions exposed to far UV radiation (Dalgarno and 
Oppenheimer 1974). It was then considered by Barsuhn & Walmsley (1977) and Clavel et 
al. (1978), who studied the chemical and thermal equilibrium in dark clouds exposed to far 
UV radiation and cosmic rays. Clavel et al. explicitly included the contribution of every 
reaction in their chemical network to heating. Chemical heating has also been widely used 
in the study of planetary atmospheres (e.g., Roble et al. 1987). 

The results of this paper are based on the fact that roughly half of the energy generated 
by cosmic rays and X-rays comes goes into ionization of the gas and roughy half into its ex- 
citation, more specifically according to the way the individual processes are treated by DYL. 
Roughly half of the gas heating comes from the reactions of various ions with the assumed 
mainly neutral gas. Potentially an equal amount of heating can arise from dissociation and 
from rotational and vibrational excitation, but the yield from excitation depends on whether 
the physical conditions are conducive to the quenching of the excited levels. If quenching is 
inefficient, the levels decay with the emission of radiation, which can escape or be absorbed 
by dust. This possibility is even more important for the excitation of the singlet levels of H 2 , 
e.g., the B and C levels, and the excitation of the n = 2 level of H in atomic regions. We 
have not attempted to follow the fluorescent radiation to determine how much escapes and 
how much is absorbed. This is an important issue when considering the broader thermal 
properties of the gas. The treatment of the fluorescent radiation involves radiation transfer 
and depends on the properties of the dust, i.e., it is application specific. In contrast, our 
goal has been to treat one relatively well-defined part of the thermal problem of interstellar 
and circumstellar molecular gas. 

In pursuing this goal, we have also ignored the direct interaction of X-rays and cosmic- 
rays with the dust, of some interest because of the possibility that they may release fast 
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electrons from the dust. Although this occurs, it is relatively unimportant. First of all, 
the dust cross section per H nucleus for the MRN distribution for the diffuse ISM is 1.6 x 
10 _21 cm~ 2 . For molecular clouds and disk atmospheres, it will be even smaller. The main 
inelastic cross sections for a keV electron with H 2 are typically several times lCT 17 cm -2 , 
so the probability that an incident X-ray or cosmic ray interacts with dust is roughly 1000 
time less than with H 2 . And when this rare event does occur, no more than 1% goes into in 
photoelectrons, as shown by Dwek and Smith (1996) for the EUV/Xray bands. 
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A. Chemical Heating in Atomic Regions 

Table A provides the same information on the energy partitions for an atomic H, He 
mixture as Table 1 does for a H 2 , He mixture. Although less energy is expended in ion 
production (because of the smaller ionization potential of H compared to that for H 2 ), the 
total number of ions is essentially unchanged at 27.2 ion pairs per keV. 



This manuscript was prepared with the A AS IATj^X macros v5.2. 
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Table A. DYL Energy Partitions for a He, H mixture, x e = and E — 1 keV 



Process 


W 


1 keV/W 


Specific energy 


Total energy 




(eV) 




(eV) 


(eV) 


TT_I_ 

H+ 


39.8 


25.1 


13.60 


342 


He + 


487 


2.05 


24.6 


51 


He 2+ 


16400 


0.061 


54.4 


3 


Total Ion production 








396 


H2 1 S 


267 


3.75 


10.2 


38 


H2 X P 


33.1 


30.2 


10.2 


308 


En = 3 


191 


5.24 


12.1 


6 


Hn>3 


236 


4.24 


13.2 


56 


Total H excitation 








408 


He2 1 S 


17300 


0.058 


20.6 


1.2 


He2 1 P 


2080 


0.481 


21.2 


10.2 


t t n3o 

He2 J S 


50700 


0.020 


19.8 


r\ on 

0.39 


rie Z, r -\-n > z, 


oUoUU 


U.UoZ 


oi n 

Zl.U 


U.Oo 


Hen > 2 1 


3790 


0.264 


23 


6.1 


Total He excitation 








19 


Total excitation 








427 


Heating 








117 


Grand total 








940 



The heating of an atomic region is simpler to treat than a molecular region because of 
the absence of rotational and vibrational excitation. The direct heating by elastic scattering 
is usually important. We ignore the radiation emitted following the electronic excitation of 
the levels of H and He in Table A, and focus here on the chemical heating. 
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In addition to radiative recombination, 

e + H + ->H + /w/, a(H + ) ~ 2.12 x 10" 10 r- - 73 cm 3 s _1 , (Al) 

H + can be destroyed by radiative attachment to H to form H^~, 

H + + H -> H+ + hv, k x& ~ 2.10 x 1(T 23 T L80 cm 3 s~ 1 . (A2) 

The rate coefficient for radiative association is especially small, and the probability to form 
this way is given by, 

P(H + , H + ; H) = fcraX ( H ) = X ( H ) _ (A3) 

fc ra x(H) + ax e x(H) + (a/k ra )x e 

In order for radiative recombination to be important, the electron fraction must approach 

the value, 

^^ 10 -13 T 2.53_ (A4) 

a 

This requirement is relatively easy to satisfy, e.g., k ra /a ~ 1CT 8 at 100 K and ~ 10~ 4 at 
4000 K. 

The energetics of the two processes are different. In both cases the radiation is lost or 
absorbed by dust grains and will be ignored. In an atomic H region the newly formed 
can charge-exchange with H to form H 2 , 

H+ + H^H 2 + H+, k = 6.4 x 10" 10 cm 3 s _1 . (A5) 

This reaction can provide up to g(H + ; H) = IP(H.2) — JP(H 2 f ) = 1.83 eV in heating, where 
IP stands for ionization potential, but it has to compete with dissociative recombination, 
Eq. [T71 which has a branching ratio, 

_ (3x e _ (f3/k)x e , . 

e ~ (3x c + fex(H) ~ (f3/k)x e + x(H) ' 1 ' 

The ratio k//3 — 5.33 x 10 _3 T ' 4 determines the electron fraction where dissociative recom- 
bination is important, i.e., in regions where x e ~ 0.01. 
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We follow Section 3 and express the heating per ion pair for atomic regions due to the 
H + ion as, 

Q chem (H+; H) = m fi H) m s F(H+; H) P(H+ H+; H) [B u g H (H+; H) + B c g c (H+; H)] , 
x(H) + x{H. 2 ) 

(A7) 

where P(H + ; H) is the fraction of all ions that are H + ; P(H + , H^; H) is the probability that 
an H + ion is converted into H^~; gn(H + ;H) = 1.83 eV is the net heating from reactions 
Eqs. IA2I and IA5( Ph is the branching ratio for reaction Eq. \A5\ 



fcs(H) fcx(H) 
H kx(E) + (3x e x(R) + (P/k)x e 1 ; 

and q e (H + ; H) = 11.0 eV is the net heating from radiative charge exchange, Eq. \A2\ followed 
by dissociative recombination, Eq. [T71 From Table A, we find that the fraction of H + ions 
P(H + ; H) = 0.92 for the parameters in that table (E — 1 keV and x c = 0); it increases slowly 
with decreasing electron energy. 

The maximum chemical heating per ion pair coming from H + depends on how the it 
generates is destroyed: 1.7 eV by charge exchange and 10.1 eV by dissociative recombination, 
with the latter only occurring for relatively large electron fractions. However in most HI 
regions, the dominant factor is the probability of forming H + , P(H + , H^; H) in Eq. IA71 which 
can greatly reduce Q c hem(H + ; H), simply because radiative association in this case is such a 
weak process. Thus chemical heating in atomic H regions is usually negligible. 

In addition to radiative recombination, 

e + He+ ->• He + hv , a(He + ) ~ 2.12 x 10- 10 T-°- 73 cm 3 s _1 , (A9) 

He + can be destroyed by radiative charge exchange, 

He + + H -> H+ + He + hv , A; chex ~ 1.6 x 10~ 16 T - 50 cm 3 s _1 , (A10) 

and by radiative association to form HeH + , 

He + + H -> HeH+ + hv , fc ra ~ 1.33 x 10" 14 r ' 37 cm 3 s _1 . (All) 
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All of these radiative processes are weak, although the last two are not as weak as the 
corresponding reactions just discussed that start from H + + He. These reactions have the 
potential to provide greater heating because of the larger energy of the He + ion, except that 
much of the available energy is lost in radiation. We therefore ignore radiative recombination 
and focus on radiative charge exchange, Eq. IA101 and radiative association, Eq. IA11I The 
former process generates IP (He) — 7P(H) = 11.0 eV in recombination radiation rather than 
gas heating. The heating from the formation of the HeH + ion, discussed in the previous 
section, has a maximum value 1.83 eV. But this value must be reduced by the fraction of 
He + which, according to Table A for E = IkeV, is _F(He + ,H) = 0.075. It may also be 
reduced by the branching ratio for radiative association, 

5 ra (He+) = , , (A12) 

"'ra ~r "-chex t OlX e 

which is often less than unity. For example, in warm regions charge exchange is more 
important than radiative association, and radiative recombination is competitive with the 
other reactions for x e > 10~ 3 . Thus we cannot expect the chemical heating from He + to be 
any more than 0.14eV, and thus negligible, as we found also for the H + ion. 



